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The main thrust of our research efforts during the period of the 
NASA grant involved theoretical studies of the electronic structure of 
metal clusters, in particular clusters of Group IIA and IIB atoms. Early 
in the project it became clear that electron correlation involving 
d orbitals plays a more important role in the binding of these clusters 
than had been previously anticipated. This necessitated that we replace 
our computer codes for calculating two electron integrals and for constructing 
the resulting Cl Hamiltonions with newer, more efficient procedures. Hence, 
considerable time was spent on program modification, interfacing, and 
testing. (See Appendix A) 

I. Group II A metal atom clusters 

Perhaps the most significant finding resulting from our studies is 

2 

that although Mg^ is essentially unbound at an SCF level of treatment the 
inclusion of electron correlation effects involving d functions gives rise 
to appreciable (^0.8 eV) binding.^ This is the first example of strong 
binding originating from correlation effects alone. This work was per- 
formed in collaboration with Professor C. Dykstra and utilized the CEPA 
and SCEPA methods for treating electron correlation. These calculations 
on Mg^ included explicitly all electrons and represents one of the largest 
(in terms of the number of electrons treated explicitly and number of 

basis functions employed) calculations on a metal cluster to date, 
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During the course of our study of Mg^, Pople and coworkers published 
a paper on the electronic structure of beryllium clusters. These authors 
also found that correlation effects involving d functions plays an 
important role in the cluster binding. However, a fundamental difference 
exists between the Be^ and Mg^ cases in that Be^ displays appreciable 
binding even in the absence of correlation effects. The discoveries that 
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d functions play a fundamental role in the binding of beryllium and 

magnesium clusters has led us to perform new calculations on Be 2 and 

Here we have found that it is essential to Include both d functions and 

correlation effects to obtain the correct ordering of the pi and sigma 
4 

anion states. 

We are presently carrying out a new set of calculations on Mg^ 
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utilizing a local potential to model the Is 2s 2p atomic cores. If 
these local potential calculations successfully reproduce the results 
of the all-electron calculations, then, we can confidently employ local 
potentials to study the larger clusters. 

As part of our continuing effort to understand electron attachment 
to group IIA metal clusters, we have also calculated^ electron affinities 
of atomic Be, Hg, and Ca. For the atoms the anion states lie energetically 
above the ground states of the neutral atoms, requiring scattering calcu- 
lations for their investigation. Using the Harris-Michels method, we 
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have shown that all three atoms have low-lying P anions. For Mg a broad 
2 

feature due to D Mg is found near 2.5 e.V in the scattering cross 

2 

section. For Ca the D anion is found to lie only 0.7 ev above the 
P anion. 

II. Cu^ and CuH . 

During the course of the grant we began calculations on Cu^ and CuH 
6 

and their anions . In these studies we have utilized the effective 

potential of Topiol et .al .^ for the copper atoms. The first challenge 
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was to find a compact Cu basis set suitable for the d s and d s states 

10 2 

of Cu as well as the d s state of Cu . After much experimentation, we 
arrived at a [4s3p5d/4s3p2d] basis set. With this basis we are able to 
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obtain a SCF -»■ d^s^ excitation energy of 0.59 eV in fair 

agreement with the all-electron value of 0*34 eV and a SDCI (single and 

double excitation configuration interaction) excitation energy of 1.41 eV 

in excellent agreement with the experimental value of 1.49 eV. A SDCI 

with this basis yields a dissociation energy of Cu 2 of 1.87 eV compared 

with the 2.05 eV experimental value. The above results are noteworthy 

since other investigators have concluded that a 5d/3d contraction is 

9 2 10 

required to describe the d s ^ d s atomic excitation energy and pre- 
sumably also to describe the bonding in the CU 2 molecule. The calculations 
on the CU 2 and the CU 2 potential energy curves as well as on CuH and 
CuH are near completion. We expect to submit this work in the near future 
to the Journal of Chemical Physics for publication. 

III. Reactions of Metal Atoms with H ^Q* 

The second major thrust of our research program has been the investi- 
gation of the reactions of metal atoms with small molecules, especially 

H 2 O. In particular, we have characterized the M*0H2 and HMOH species 
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for M= Be, B, C, Mg, Al, Si. * For each of these species the geometries, 
stabilities, and vibrational frequencies have been determined. Perhaps 
the most interesting finding of these calculations is that for the B, C, 

Al, and Si insertion products both cis and trans isomers exist. Experi- 
mental evidence for the existence of cis/trans isomers of HSiOH has 
recently been found . 

IV. Other Work Carried Out Under the NASA Grant 

Theoretical investigations of several other systems were also com- 
pleted during the course of our NASA grant. Using 2^^ order many-body 
perturbation theory we have determined the polarizabilities of LiF and 
BeO, BF, and LiCl. The results for LiF and BeO and those for 
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BF have been published. Those for LiCl are near completion. The 

calculations on LiCl are of particular interest since in this case we 

have determined several properties using different pseudopotentials for 

the chlorine atom. Finally, we have utilized^^ both the configuration 

interaction and many-body perturbation theory (through fourth-order) 

to calculate the electron affinities of and C^. We plan to extend 

these studies to include larger carbon clusters. 
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Appendix A : Program Implementation and Development 

During the course of the NASA grant we Implemented the Pltzer-Schaefer 
Integral program, which Is considerably faster for Integrals over d functions 
than the MOLECULE program we had been using. A significant gain In speed 
resulted from a switch over to Shepard's Unitary-group multlconflguratlonal 
self-consistent field (MCSCF) and configuration Interaction (Cl) programs. 
These programs were modified to allow for the use of local potentials to 
model the core electrons and to use Input vectors from Purvis' SCF program. 
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